home *** CD-ROM | disk | FTP | other *** search
/ Cream of the Crop 26 / Cream of the Crop 26.iso / os2 / octa209b.zip / octave-2.09 / SCRIPTS.ZIP / scripts / linear-algebra / trace.m < prev    next >
Text File  |  1996-07-12  |  1KB  |  40 lines
  1. ## Copyright (C) 1996 John W. Eaton
  2. ##
  3. ## This file is part of Octave.
  4. ##
  5. ## Octave is free software; you can redistribute it and/or modify it
  6. ## under the terms of the GNU General Public License as published by
  7. ## the Free Software Foundation; either version 2, or (at your option)
  8. ## any later version.
  9. ##
  10. ## Octave is distributed in the hope that it will be useful, but
  11. ## WITHOUT ANY WARRANTY; without even the implied warranty of
  12. ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13. ## General Public License for more details.
  14. ##
  15. ## You should have received a copy of the GNU General Public License
  16. ## along with Octave; see the file COPYING.  If not, write to the Free
  17. ## Software Foundation, 59 Temple Place - Suite 330, Boston, MA
  18. ## 02111-1307, USA.
  19.  
  20. ## usage: trace (x)
  21. ##
  22. ## Returns the trace (the sum of the diagonal elements) of x.
  23.  
  24. ## Author: jwe
  25.  
  26. function y = trace (x)
  27.  
  28.   if (nargin != 1)
  29.     usage ("trace (x)");
  30.   endif
  31.  
  32.   [nr, nc] = size (x);
  33.   if (nr == 1 || nc == 1)
  34.     y = x(1);
  35.   else
  36.     y = sum (diag (x));
  37.   endif
  38.  
  39. endfunction
  40.